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We present a new normalization condition for the axial current, which is derived from the PC AC relation 
with non-vanishing mass. Using this condition reduces the O(rom) corrections to the axial current normalization 
constant Za for an easier chiral extrapolation in the cases, where simulations at zero quark-mass are not possible. 
The method described here also serves as a preparation for a determination of Zk in the full two-flavor theory. 
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1. Introduction 

Due to the explicit breaking of chiral symme- 
try in a lattice theory with Wilson quarks a finite 
renormalization constant Zp^ for the isovector ax- 
ial current is needed so that Ward identities take 
their canonical form up to small lattice correc- 
tions. Otherwise a safe extrapolation to the con- 
tinuum can not be achieved for physical quanti- 
ties involving this current. 

The function ZpXso) Q currently used by the 
ALPHA collaboration was obtained by requiring 
that certain chiral Ward identities are satisfied in 
the 0(a) on-shell improved lattice theory up to 
O(a^) cutoff effects. This normalization condition 
is evaluated at zero quark mass and hence the 
mass term in the PCAC relation was dropped. 

However, even in the framework of the 
Schrodinger functional simulations at vanishing 
quark mass are not possible if the physical vol- 
ume, or equivalently the renormalized coupling, 
is too large. For the large couplings one there- 
fore evaluates the normalization condition at fi- 
nite quark mass and extrapolates to the chiral 
point. Since the mass term was ignored in the 
derivation of the normalization condition this in- 
troduces an error of Oir^rri) in Zp^. The Sommer 
scale To is used as a typical hadronic scale. 

Here we present a normalization condition for 
the axial current, which is derived using the full 
PCAC relation. This reduces the O(rom) errors 
in Za to Oiara) and therefore results in a much 
flatter chiral extrapolation. 



2. The axial Ward identity 

By performing local infinitesimal symmetry 
transformations of the quark and anti-quark 
fields in the Euclidean functional integral one de- 
rives the Ward identities associated with the chi- 
ral symmetry of the continuum action. An ax- 
ial transformation of the isospin doublet quark 
fields gives the partially conserved axial current 
(PCAC) relation 



kd^Al(^^)0) - 2m{P\x)0) = -{SO), 



(1) 



where SO is the variation of some operator O. 
A'^{x) and P"(x) are the isovector axial current 
and isovector pseudoscalar density, respectively. 

Let TZ be the space-time region where the chi- 
ral transformation is applied. The axial current 
Al, (y) is inserted as an internal operator (y 6 7?,) . 
If Ocxt denotes a polynomial in the basic fields 
outside this region the integrated form of the ax- 
ial current Ward identity is 
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-2m y"d*a;(Jp"(x)A^(j/)axt) 

= i6'"'"(K'(2/)axt). (2) 

The right-hand side of originates from the 
variation of the internal operator A^^{y) since un- 
der chiral transformation the axial current rotates 
into the vector current V^{x). 
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3. Normalization conditions on the lattice 

A normalization condition for tlie axial current 
on the lattice is derived by demanding that eq. Q 
in terms of the renormalized currents is still valid 
up to lattice artifacts for a given choice of TZ and 

Cext- 

We use Schrodinger functional boundary condi- 
tions in 0(a)-improved lattice QCD '3l4j. i.e. we 
have an L'^ X T space-time cylinder with periodic 
boundary conditions in the spatial directions and 
fixed (Dirichlet) boundary condition in the tem- 
poral direction. 

We choose TZ =^ x [T/S, 2T/3\ and an Ooxt, 
which creates pseudo-scalar states at the upper 
and lower temporal boundary. In the renormal- 
ization scheme we employ |5I6| the normalization 
condition now gives a relation between a set of 
(improved) correlation functions on the lattice 
and the renormalization factor Za(1 + bp^amq). 

If the mass is set to zero in (j2l with our choice 
for Ooxt the lattice version of equation jSJ be- 
comes 

-^A/iA = /l) (3) 

where fi is the correlator connecting the pseudo- 
scalar boundary states and is the correlation 
function corresponding to the surface integral in 
in terms of the improved lattice currents. 
If the mass term in the PCAC relation is not 
neglected one has to evaluate an additional cor- 
relation function /p^ corresponding to the sec- 
ond term in (0). For this correlator on-shell 0(a) 
improvement breaks down due to contact terms. 
The normalization condition thus obtained is 

Zl{l + bp^am,fi^fij, - 2m/^^) = /i. (4) 

In the improved theory the normalization condi- 
tions and ^ result in an intrinsic uncertainty 
of O(a^) in Za- Condition ^ will lead to an ad- 
ditional error of O(rom) when evaluated at finite 
quark mass whereas (0} will lead to 0(am) errors 
from the contact terms in the correlator multiply- 
ing the quark mass and also from 6a, which we 
only know perturbatively. 



4. Testing the new normalization condition 

We choose the Schrodinger functional setup 
with T = |L, = and vanishing background 
field {C — C — 0) and evaluate the normaliza- 
tion conditions on 8'^ x 18 quenched lattices with 
several different values of the hopping parameter. 

Figurenshows the results from both normaliza- 
tion conditions as a function of the PCAC mass 
for (3 = 8.0 on the same set of lattices. In this 
case a simulation at the chiral point is possible 
in the Schrodinger functional. As expected the 
quark mass dependence is strongly reduced when 
the new normalization condition is used. Also 
the statistical error in Za is much smaller with 
the massive condition. 
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Figure 1. Comparison of old (massless) and new 
(massive) normalization condition at small cou- 
pling. The dotted line is the result for Zj^. 

More interesting is the case oi (3 = 6.0, where 
Za has to be extracted from a chiral extrapola- 
tion using a quadratic fit. Figure El shows that 
in this case the curvature is more pronounced 
but the mass dependence for the new method 
is still quite flat. The error in the extrapolated 
value for Za is however larger for the massive 
normalization condition: 

massive Za — 0.8135(58) , . 

massless Za = 0.8237(37) ^ ' 



3 




-0 0005 0.0000 0.0005 



0.3 r I I I 1 

0.000 0.010 0.020 0.030 0.040 

am 

Figure 2. Both normalization conditions evalu- 
ated a.t P — 6.0. The extrapolated values are 
shown in the insert. 

5. Disconnected quark diagrams 

On closer inspection the reason for this larger 
error can be found by looking at the possible Wick 
contractions for the correlation functions and 
/p^. The structure of the external operator used 
leads to disconnected quark diagrams where the 
two current insertions each couple to one tempo- 
ral end of the lattice only. 

Those cause large statistical fluctuations espe- 
cially in the volume correlator /p^, which make 
the extraction of Za from the new normalization 
condition statistically more difficult. 

These diagrams can be suppressed by using a 
different external operator in the derivation of the 
normalization conditions. This operator requires 
the existence of a third quark species (spectator) 
and can therefore be used only as an approxima- 
tion in a theory with Nf < 3. We are working on 
an extension of this argument. 

Figure 13 shows the result from both conditions 
if only connected quark diagrams are used. Note 
that in this case a linear extrapolation in the 
quark mass is possible. The extrapolated values 
for Za in this case are 

massive Za = 0.8135(20) 
massless Za = 0.8148(25). 
Here the new condition not only results in a 
flat extrapolation, it also gives an extrapolated 



value of Za, which is consistent with © and has 
a smaller statistical error. 
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Figure 3. Both normalization conditions using 
only the connected diagrams. 

6. Conclusions 

The new normalization condition including the 
mass term has only a very small dependence on 
the quark mass. If the spectator quark argument 
can be extended to the two flavor case it is there- 
fore a suitable tool for calculating the axial cur- 
rent normalization constant in the cases were sim- 
ulations at zero quark mass are not possible. 

This work is also a preparation for an un- 
quenched determination of Za- Due to algorith- 
mic problems at large values of the bare gauge 
coupling and small quark masses we will benefit 
greatly from a reliable chiral extrapolation. 
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